Muon anomalous magnetic moment 
in the standard model with two Higgs doublets 



Yue-Liang Wu and Yu-Feng Zhou 
Institute of Theoretical Physics, Chinese Academy of Sciences, 
Beijing 100080, China 

Abstract 

The muon anomalous magnetic moment is investigated in the standard 
model with two Higgs doublets (S2HDM) motivated from spontaneous 
CP violation. Thus all the effective Yukawa couplings become complex. 
As a consequence of the non-zero phase in the couplings, the one loop 
contribution from the neutral scalar bosons could be positive and nega- 
tive relying on the CP phases. The interference between one and two loop 
diagrams can be constructive in a large parameter space of CP-phases. 
This will result in a significant contribution to muon anomalous mag- 
netic moment even in the flavor conserving process with a heavy neutral 
scalar boson (m^ ~ 200 GeV) once the effective muon Yukawa coupling 
is large ( ~ 50). In general, the one loop contributions from lepton 
flavor changing scalar interactions become more important. In particu- 
lar, when all contributions are positive in a reasonable parameter space 
of CP phases, the recently reported 2.6 sigma experiment vs. theory 
deviation can be easily explained even for a heavy scalar boson with a 
relative small Yukawa coupling in the S2HDM. 
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Recently, the E821 Collaboration at BNL has reported their improved result on the 
measurement of muon anomalous magnetic moment = (g^ — 2)/2 Q. The difference 
between the measurement and the Standard Model (SM) prediction [J2| is reported 

Aa^ = a e * v - al M = 426 ± 165 x 10" 11 , (1) 

which shows a relative large deviation (2.6 a) from the SM calculation. At 95% confidence 
level, Aa M lie in the range 

113 x 1CT 11 < Aa M < 749 x 10 -11 . (2) 

As muon is about 210 times heavier than electron, it is expected that the new physics effects 
on muon anomalous moment may be considerable. A large amount of works have been made 
on checking the new physics contributions to Aa M from various models. Among those, the 
contributions from supersymmetric particles and scalar bosons as well as extra dimensions 
seem to be the most attractive ones • In this paper, we shall make a systematic analysis 
in the standard model with two Higgs doublets (S2HDM) motivated from the study of origin 



and mechanism of CP violation [^-pUj] . In this kind of model, the Higgs sector of SM is simply 
extended by including an additional Higgs doublet with all the Yukawa couplings being real. 
After spontaneous symmetry breaking, the origin of both fermion masses and CP-violating 
phases can be attributed to the well known Higgs mechanism with a single CP-phase between 
two vacuum expectation values. Namely CP symmetry is broken spontaneously [|ll|]. Thus 
all the effective Yukawa couplings in the physical basis become complex. The S2HDM 
has also been investigated from other considerations Hl2|— pTTfj . It will be seen that though 
S2HDM is the simplest extension of SM, it can provide a possible good explanation for 
the recently reported 2.6 sigma experiment vs. theory deviation of the muon anomalous 
magnetic moment. 

There are several versions of 2HDM. To avoid flavor changing neutral current (FCNC) 
at tree level, some kind of discrete symmetries are often imposed on the Higgs potential, 



this results in two kinds of models which are usually called the 2HDM of type I and II [|I8"| 
Note that by imposing the discrete symmetry, the spontaneous CP violation is forbidden 
with two Higgs doublets and all the Yukawa couplings have to be real. 

Contributions to muon g — 2 from the 2HDM of type II have been discussed in Refs. [|J . 
It seems that the one loop contribution is not large enough to explain the data. Even for 
a large value of tan/? ~ 50, one still needs a very light mass of the scalar ~ 5 GeV. If 
both the scalar and pseudo-scalar are included, due to the cancelation between them, the 
situation will be even worse. One then needs to consider two loop contributions. It was 
found H^] that from the two-loop diagrams of Barr-Zee type ||20|| , the contributions from 
pseudo-scalar is positive and could be larger than its negative one loop contributions. Thus 
the net effects for pseudo-scalar exchange become positive. Provided a sufficient large value 
of the coupling tan /3, its contributions can reach the experimental bound. By including the 
two loop diagrams, the pseudo-scalar mass must be below 80 GeV even when tan (3 is large 
and around 50 ||. To avoid the cancelation between scalar and pseudo-scalar exchange, the 
mass of the scalar boson has to be pushed to be very heavy ( typically greater than 500 
GeV). 

An alternative way of suppressing the FCNC at tree level is to take the smallness of the 
off-diagonal elements in the Yukawa coupling matrices by considering approximate flavor or 
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family symmetries [ 13 , l^ ,p-^0 . It has been shown |8|-|T0f that through the spontaneous CP 
violation, rich sources of CP violation including the KM-phase in SM can be induced from 
a single relative phase between the vacuum expectation values of the two Higgs doublets. 
In the S2HDM, as there are flavor changing scalar interactions for both neutral and charged 
scalar bosons, the bounds on the neutral Higgs mass can be released through the inclusion 
of the internal r loop at one loop level |7j. Besides this, as it will be shown below that 
the complex and flavor dependent Yukawa couplings in S2HDM may completely change the 
interference between one and two loop diagram contributions. As a consequence, one and 
two loop contributions can be both positive and provide significant contributions to Aa M . 
In the S2HDM, after spontaneous symmetry breaking, the two Higgs fields are given by 



-±(vcos/3e i5 + 0?) J' ^( us i n/ 3 + 0O) 



(3) 



where the phase 5 is the relative phase between the vacuum expectation values of the two 
Higgs doublets. This phase is the physical one as it can not be rotated away by a redefinition 
of the fermion phase. Thus in the physical basis (i.e., mass eigenstate basis), all the effective 
Yukawa couplings become complex, which leads to rich sources of CP violation. It is more 
natural to use another basis for Higgs fields through the recombination 

E x = 0i sin (3e~ tS - <p 2 cos (3, H 2 = <fi 2 cos j3e~ lS + <p 2 sin (3 (4) 

such that 

Hl= {#fl+u))' H > = (> + tf°))' (5) 

where H + is the charged scalar and R, I and H° are the three neutral scalar bosons. Here 
the neutral Higgs boson H° plays the role as the one in SM with one Higgs doublet. The 
additional two neutral scalar bosons R and / will lead to new interactions beyond the SM. 

The Yukawa interactions between scalar bosons and leptons in the S2HDM can be written 
as the sum of two parts C\ and C 2 , where L\ contains only flavor conserving interactions 
and C 2 contains flavor changing ones, i.e. 

£ Y = (V2G F ) 1 / 2 (£ 1 + £ 2 ) (6) 

with 

3 3 3 

A = V2H+ £ i h m h vir R + h°y: m h T l A + (R + H) £ ^m^A + h.c, (7) 

i i i 

C 2 = {R + l I)jzi l J L l J R + h.c, (8) 

In the above expression, ^. and stand for the flavor conserving and flavor changing 
Yukawa couplings which are in general all complex numbers. For simplicity, we will neglect, 
in the following considerations, the mixing among the three neutral scalar bosons, but their 
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masses are considered to be largely split and the mass eigenstates are denoted by h and A 
respectively. 

The phase effects of the Yukawa coupling has been discussed in the literatures 
f8|- |lT||21~lj23| , The complex couplings may lead to sizable CP violation in charged and neu- 
tral Higgs boson exchanging processes in quark sector, such as b — > s'j |9|,f2^,[22|] , K° — K°, 
B° — B°, B° s — B® and D° — D° mixings ||24j| , etc. In the lepton sector, it may result in large 
electric dipole moment of leptons pT| , pi]|J2B| . 

In the case of effective real Yukawa couplings for neutral scalar bosons, h is purely a 
scalar while A is a pseudo-scalar. The contributions to muon anomalous moment from these 
two particles always have different signs at both one loop and two loop level. Namely, for 
scalar /i-exchange, the contribution is positive from one loop and negative from two loop. 
But for pseudo-scalar A-exchange , the one loop contribution is always negative and the two 
loop contribution is positive. As the experimental data indicate that the deviation of g — 2 
from SM is positive. The large two loop contribution may lead to the conclusion that the 
A-exchange must be dominant and the /i-exchange must be made to be negligible small. 

However, in the general case, the situation can be quite different. As all the couplings 
are complex in S2HDM, the Yukawa couplings for neutral scalar bosons h and A contain 
both scalar and pseudo-scalar type fermion interactions. Further more, the couplings are all 
flavor dependent, which is far away from the 2HDM of type I and I, where the couplings 
are given by a simple function of tan (3. The fact that one loop and two loop contributions 
depend on different complex Yukawa couplings provides an alternative possibility that they 
may not always cancel each other. In a large parameter space, the two kind of contributions 
may be constructive and result in a large value of muon g — 2 

The above discussion is also suitable for SUSY model [|5|,|6[] as the Higgs sector 
of it is exactly like the 2HDM of type II. Recently, the contribution from two 
loop diagrams with charged Higgs in SUSY is also discussed ||27|| . However, In 
the case of SUSY, the main contribution may arise from the chargino-sneutrino 



and neutralino-slepton couplings at one loops |[28|| . 

The one loop flavor conserving contributions from the scalar and pseudo-scalar have been 
investigated in Refs. pOf. In S2HDM, the results are given by 
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where g is the weak gauge coupling constant. In the limit of m^/m^ -C 1, the one loop 
results can be approximately written as 
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where 5^ is the phase of £ M with ^ = |£ M |e 4<5 '\ The one loop flavor conserving contribution 
to Aa M is plotted in FigPJa. It can be seen that its contributions decrease rapidly with m^ 
increasing, which means that the one loop diagram cannot give large contributions to Aa M 
except for a very small value of m^ ~ 10 GeV or a very large value of ~ 70. Note that 
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in the S2HDM, the contributions from h and A can be negative and positive depending on 
the sign of cos25 M , which is completely different from the other type of 2HDM, such as type 
I and type E, where the h(A) loop diagrams always give positive (negative) contributions. 

In the case of one loop flavor changing processes ( see Fig.[l] b), the loops with internal 
r-lepton play an important role as it is much heavier than the /i-lepton. This then leads to 
an enhancement factor of the order m 2 /m 2 l ~ (9(10 2 ) relative to the flavor conserving one. 
In S2HDM, the contributions from h(A) exchange are given by 

AoJW = ±4-i^- f'"% - f] M 2 cos2V (12) 
" 32tt 2 m 2 w m 2 h(A) \ m 2 2 J 

In obtaining the above expression, we have taken £ Mr = £ rM for simplicity. The flavor 
changing contributions to Aa M are presented in Fig.|2|b **with different values of |£ MT |. The 
figure shows that the contributions to Aa M may be considerable large when |£ Mr | is large. 
For |^ T | = 30 ~ 50, the recently reported 2.6 sigma experiment vs. theory deviation can 
be easily explained even for a heavy scalar boson rrih > 100 GeV. So far, there are no strict 
constraints on the values of the coupling |£ MT |. Studies on the rare decays r — ■> //(e)7 and 
t — > 3/i(e) and the electric dipole moment of r will be useful to provide an interesting 
constraint on the parameter. However, as the relevant experimental data at present are 
primitive and such processes often contain more couplings such as £ r and £ Te , the resulting 
constraints can not be clearly obtained and they are not yet be very strong. To obtain the 
upper bounds on £ Tfl , further studies are needed. 

Now let us discuss the two loop effects on Aa M . Naively speaking, relative to one 
loop diagram, the two loop diagram will receive an additional suppression factor 
of a/An ~(9(10 3 ) and is in general small. However, there exists an special kind 
of two loop diagrams which can give significant contributions through the large 
Yukawa couplings between Higgs bosons and heavy fermions. These processes 
which are often referred as Barr-Zee mechanism are shown in Fig.(|l].)c and d . 
In this kind of diagrams, the Yukawa coupling will contribute a enhance factor 
of rrif/m^. if / is heavy quark such t— quark, it can reach the order of O(10 3 ). 
Thus the two loop contribution could be sizeable. The large Yukawa couplings 
of order O(50) is need for explain the mass ratio between top and bottom quark 
in GUT ||29|| . Note that arbitrary large Yukawa couplings may lead the theory 
to be non-perturbative . This requires that the couplings have to be less than 
l/(v / 2(v / 2C/) 1 / 2 m/), where Gj is the Fermi constant. For example, for t(6)-quark, 
it is less than 1(44). In the virtual loop corrections, as there are additional loop 
suppression factor of (47r) 2 . The constraints may be weak. Further more, the 
loop integrals always decrease with Higgs mass growing, and the Higgs loop 
contribution will decouple at large Higgs mass limit. 

For the case that the fermion loop is dominated by top quark, the Barr-Zee type two 
loop diagrams may lead to a sizable Aa^. In S2HDM, their contributions are given by 



2 ^,2 lyYl ^ 
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where iV c = 3 is the color number and qt = 2/3 is the charge of top quark, a = 1/137 and 
sin 2 #^ — 0.23 are the fine structure constant and weak mixing angle respectively. f(z) and 
g(z) are two integral functions which have the following form py] 

/(*) = -z t d X l ~ 2x{ \- x) In ^-^ (15) 
J v ' 2 Jo x(l - x) - z z y ' 

t s 1 f 1 , 1 , x(l—x) 

g(z) = -z dx— r ln-^ '- 16 

yv ; 2 Jo x(l-x) -z z K J 

for large z, f(z) ~ 1/3 In z + 13/8, g{z) ~ 1/2 In* + 1 and /(l) ~ 1/2, #(1) ~ 1. The 
pure two loop contributions are plotted in Fig.|2|b. Unlike in the one loop case, where Aa M 
decreases rapidly with growing rrih, the two loop contributions decrease relative slowly and 
their signs depend on the value of 5^ T . Therefore for a very large value of rrih, the two 
loop effects become dominant. Another difference is that in the one loop case the new 
physics contributions only depend on while the two loop contributions depend on two 
couplings £ M and £ t , if the fermion loops for b or r are included |||6||, it will depend on more 
parameters. As the two couplings are complex numbers, the interference between one and 
two loop diagrams may not always be destructive. There exists a large parameter space in 
which the one and two loop contributions are all positive. Thus they can result in a large 
contribution to Aa M . The constraint of |£t| has been studied in B° — B° mixing and b — > s^y 
as well as the neutron electric dipole moment [p^,3,|10|, pl] , p5|1 , the typical absolute value for 



|£t| is of the order 0(1). Taking \£ t \ = 1 and 5 t = as an example, the sum of one and two 
loop contributions from h— scalar exchange reads 



2 



a l£J m u. 

Aa = — — 

11 sin 2 6w nTw 



(i7) 



It is clear that if 5^ lies in the range 37r/4 < 5^ < 5n/A, the one and two loop contributions 
will be constructive. The detailed numerical results for two different values of 8 t are plotted 
in Fig.|3| and Fig.||, where 5^ runs in a wide range. It can be seen that for large values of 
5^ ~ 7r, the contributions from h— scalar exchange can reach the experimental lower bound 
even for a heavy scalar with rrih ~ 200 GeV. This is quite different from the existed 2HDM 
calculations in the literature, where the allowed range for the scalar boson A must be less 
than 100 GeV , and the scalar boson h must be much heavier than A, so that its negative 
contributions are negligible. 

In a large parameter space, it is also possible that contributions from one loop flavor 
conserving and one loop flavor changing as well as two loop diagrams are all positive. In 
this case, the current data can be easily explained even with a relative small value of |^ r |. 
The results for such a situation are plotted in Fig.[5[ It is seen that in this case a reasonable 
value of |^ T | ~ 10 is large enough to reach the experiment bound. 

In conclusion, we have investigated all possible new contributions to the muon anomalous 
magnetic moment from scalar interactions in S2HDM. Though S2HDM is regarded as the 
simplest extension of SM, it can provide a possible good explanation on the recently reported 
2.6 sigma experiment vs. theory deviation. In particular, the important effects of the CP 
phases appearing in the Yukawa couplings have been studied. Unlike the 2HDM of type I and 
E, the two loop h— scalar contribution to Aa^ could be positive. The interference between 
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one and two loop contributions can be constructive and result in a considerable contribution 
to Aa M , so as to coincide with the recent E821 experiment data. In the case of constructive 
interference, the experimental data Aa M can be understood even for a heavy scalar boson 
h with mass m h ~ 200 GeV, once the Yukawa coupling is large, ~ 50. In general, 
contributions from one loop flavor changing diagrams are likely to be the dominant one in a 
large parameter space. Especially, in the case that contributions from all diagrams, namely, 
the one loop flavor conserving and one loop flavor changing as well as two loop diagrams 
in S2HDM with a reasonable parameter space, are all positive, one can easily explain the 
current experimental data of the excess of the muon anomalous magnetic moment. 
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FIGURES 
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FIG. 1. Feynman diagrams contributing to Aa M in S2HDM. 

(a) one loop flavor conserving. 

(b) one loop flavor changing. 

(c) and d) two loop diagrams of Barr-Zee type 
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FIG. 2. 

a) Single one loop flavor conserving contribution to Aa M as function of m/,. The solid, dashed and 
dotted cures correspond to 1^1=30,50,70 respectively. The two dot-dashed horizontal lines present 
the current experimental bound at 95% confidence level 

b) Single two loop contribution with being fixed at 50 and 5^ varies from 0, 7r/3, 2-7r/3 to 7r( 
from down to up). 
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FIG. 3. The sum of the one loop flavor conserving and the two loop contribution to Aa^ as 
function of m^. for \^ t \ = 1 and St = 0. 

(a) The solid, dashed and dotted cures correspond to =30,50,70 respectively, with 5^ = 0. 

(b) The same as (a) with <5 M = tt/3. Note that all the contributios are negative. 

(c) The same as (a) with <5„ = 2tt/2>. 

(d) The same as (a) with 5^ = it. 
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FIG. 5. 

(a) The single flavor changing contributon, the solid, dashed and dotted curves corresponds to 
|^ T | = 10, 30,50, with V = 0. 

(b) The sum of one loop flavor conserving, one loop flavor changing and the two loop contributions 
with |£ t | = 1, St = 7r/3 and = 30, <5 M = 2ir/3. The dot-dashed, solid, dotted and dashed curves 
correspond to |^ r | = 5, 10, 20, 30, with 5^ T = 0. 
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